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Abstract 

Background: Microbe-associated molecular patterns, such as those present in bacterial flagellin, are powerful 
inducers of the innate immune response in plants. Successful pathogens deliver virulence proteins, termed 
effectors, into the plant cell where they can interfere with the immune response and promote disease. Engineering 
the plant immune system to enhance disease resistance requires a thorough understanding of its components. 

Results: We describe a high-throughput screen, using RNA sequencing and virus-induced gene silencing, to identify 
tomato genes whose expression is enhanced by the flagellin microbe-associated molecular pattern flgll-28, but 
reduced by activities of the Pseudomonas syringae pv. tomato (Pst) type III effectors AvrPto and AvrPtoB. Gene 
ontology terms for this category of Flagellin-induced repressed by effectors {FIRE) genes showed enrichment for genes 
encoding certain subfamilies of protein kinases and transcription factors. At least 25 of the FIRE genes have been 
implicated previously in plant immunity. Of the 92 protein kinase-encoding FIRE genes, 33 were subjected to 
virus-induced gene silencing and their involvement in pattern-triggered immunity was tested with a leaf-based 
assay. Silencing of one FIRE gene, which encodes the cell wall-associated kinase SIWAKl, compromised the plant 
immune response resulting in increased growth of Pst and enhanced disease symptoms. 

Conclusions: Our transcriptomic approach identifies FIRE genes that represent a pathogen-defined core set of 
immune-related genes. The analysis of this set of candidate genes led to the discovery of a cell wall-associated 
kinase that participates in plant defense. The FIRE genes will be useful for further elucidation of the plant immune 
system. 



Background 

The plant immune system involves two related inducible 
responses. The first response is activated by the detec- 
tion of microbe-associated molecular patterns (MAMPs) 
by the extracellular domains of pattern recognition re- 
ceptors (PRRs) [1]. A small number of PRRs have been 
identified in plants that recognize MAMPs derived from 
flagellin, elongation factor Tu, an ethylene-inducing 
xylanase, and certain non-proteinaceous MAMPs [2]. 
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Additionally, damage-associated molecular patterns 
(DAMPs), which typically appear in the apoplast as a 
consequence of pathogen attack, function as host-derived 
elicitors [3]. The activation of pattern-triggered immunity 
(PTI) by MAMPs and DAMPs leads to changes in the 
intracellular calcium concentration, production of reactive 
oxygen species, activation of mitogen-activated protein kin- 
ase (MAPK) cascades and transcriptional reprogramming 
[4]. These events lead, in a largely unknown manner, to 
inhibition of pathogen growth and suppression of disease. 
Successful pathogens deliver virulence proteins (effectors) 
into the plant cell and a majority of these proteins appear 
to function by interfering with host immunity-associated 
events triggered by MAMP recognition [5]. A second plant 
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defense response, effector-triggered immunity (ETI), can be 
activated in cases where a specific pathogen effector or its 
activity is recognized by a host nucleotide-binding leucine- 
rich repeat-containing (NB-LRR) resistance protein [1]. 

The interaction of tomato with Pseudomonas syringae 
pv. tomato (Pst) is a powerful model system for under- 
standing the molecular basis of bacterial pathogenesis 
and plant immunity [6]. Pst enters the apoplastic space 
of tomato leaves through stomata or wound sites. The 
plant responds to Pst even at the stomatal entry stage, 
although the bacterium can overcome this response with 
the aid of the phytotoxin coronatine [7]. Pst present in 
the apoplast is detected by PRR-mediated recognition of 
various MAMPs thereby activating PTI. The best char- 
acterized of these MAMPs is a 22-amino acid region of 
flagellin (flg22), which is recognized by FLS2, a PRR with 
extracellular leucine-rich repeats (LRRs) and an intracellu- 
lar protein kinase domain [8]. The mechanism by which 
FLS2 activates intracellular events involves another LRR 
receptor-like kinase, BAKl, as well as cytoplasmic protein 
kinases of the PBSl-like family [9-11]. Recently, another 
MAMP derived from flagellin, flgII-28, has been identified 
[12]. The PRR that detects flgII-28 is unknown, but it 
appears to play an important role in activating PTI in 
solanaceous species [13]. 

Two Pst effector proteins, AvrPto and AvrPtoB, have 
been studied extensively and found to play multiple im- 
portant roles during the tomato-P5^ interaction [14]. 
These roles are partially redundant and deletion of both 
effector genes is required before a marked decrease in 
virulence is observed. Some of this redundancy is attrib- 
utable to the fact that each effector binds and interferes 
with protein kinase domains of the FLS2-BAK1 com- 
plex, thereby disrupting the host response to flg22 [14]. 
However, each effector targets additional host proteins 
using independent domains found in each effector 
[15,16]. In a weU-studied example of ETI, the host 
protein kinases Fen and Pto appear to act as decoys 
of the real kinase virulence targets. Interaction of 
AvrPto and AvrPtoB with these decoys triggers a host 
immune response mediated through the NB-LRR pro- 
tein Prf [17]. 

Host responses associated with PTI and ETI are complex 
and have been studied by both reductionist approaches 
focused on individual components and mechanisms as 
well as by large-scale systems biology approaches look- 
ing at the totality of metabolomic, proteomic and tran- 
scriptomic changes. Previous transcriptomic studies in 
Arabidopsis thaliana, tomato and other species have 
relied primarily on microarrays and cDNA-amplified 
fragment length polymorphism (AFLP) methods, which, 
though limited in their dynamic range and, in some cases, 
in their genome coverage, have yielded important in- 
sights into PTI, ETI and the role of pathogen effectors 



in suppressing the former responses [18-25]. However, the 
emergence of RNA sequencing (RNA-seq) and associated 
statistical analyses as well as the availability of improved 
gene annotation now allows for the development of a 
much more comprehensive view of subtle gene expres- 
sion changes, and how pathogen effectors impact these 
changes as part of their immunity-suppressive activities. 
Here we have taken advantage of the recently completed 
tomato genome sequence [26], the comprehensive nature 
of RNA-seq analysis, virus-induced gene silencing (VIGS) 
and the well- characterized tomato-P5^ system to develop 
a genome-wide screen for genes whose expression is spe- 
cifically altered in response to both MAMP detection and 
subsequent activities of pathogen effectors. 

Results 

Development of an RNA sequencing approach to 
examine the plant response to MAMPs present in various 
bacterial pathogens and non-pathogens 

To characterize the PTI response of tomato and the ef- 
fect of the delivery of a subset of effectors, we performed 
an RNA-seq analysis of tomato Rio Grande prf3 leaves 
challenged with either the flgII-28 peptide or the bacter- 
ial strains described in Table 1. We focused on flgII-28 
because it has recently emerged as a second MAMP 
derived from flagellin, in addition to flg22, which is 
perceived specifically by solanaceous species [12,13]. 
In this case we chose 1 \\M concentration, which is the 
same as the lowest amount used for flg22 treatments in 
previous microarray studies oi Arabidopsis [19,23]. Using 
Pseudomonas fluorescens, Pseudomonas putida and flgll- 
28 we investigated PTI-associated transcriptional changes, 
while Agrobacterium tumefaciens and Pseudomonas 
syringae pv. tomato {Pst) strain DC3000 AhrcQ-UAfliC 
were chosen to examine the effects due to MAMPs other 
than flagellin. The comparison between the response 
to DC3000 and the mutant DC3000 AavrPtoAavrPtoB 
was designed to uncover the effect of two 'early-acting' 
type III effectors [27]. We chose a sampling time (6 h) 
and a bacterial titer (5 x 10^ cfu/mL) to capture early 
PTI-associated transcriptional changes and also the effects 
of early-acting effectors on these changes. For both the 
non-pathogenic Pseudomonas species and A. tumefaciens 
we used the same time but a higher titer (10^ cfu/mL) to 
maximize the plant response. 

The RNA-seq analysis over all treatments revealed 
reads per Idlobase of exon model per milUon mapped 
reads (RPKM) ranging from 0 to approximately 39,000 
for the 34,727 genes predicted in tomato [26]. We focused 
on genes with >3 RPKM in at least one treatment based 
on analysis of the RPKM of genes known to play a role 
in immunity (Additional file 1: Figure SI). We refer to 
this group of 17,597 genes as the expressed genes. All 
of our RNA-seq data are available from the Tomato 
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Table 1 Details of the treatments used for RNA sequencing analysis of the tomato immune response 



Treatment 


Concentration 


Comment 


flgll-28 


1 mM 


28-amino acid peptide from Pseudomonas syringae pv. tomato Tl 


P^pi irlnmnnn^ wrinnnp nv tnmntn (P^f) DC^OOC) 


5x10^ cfu/mL 


Wild-tvne P<^t <^train 


DC3000AhrcO-UAfliC 


5x10^ cfu/mL 


DC3000 mutant, lacl<ing a functional T3SS and flagellin 


DC3000AavrPtoAovrPtoB 


5x10^ cfu/mL 


DC3000 mutant, witli a functional T3SS, but lacking avrPto and avrPtoB 


Pseudomonas fluorescens 55 


10^ cfu/mL 


Soil bacterium 


Pseudomonas putida 1^2240 


10^ cfu/mL 


Soil bacterium 


Agrobacterium tumefociens GV2260 


10^ cfu/mL 


Disarmed plant pathogen 


MgCb 


10 mM 


Mock treatment 






Untreated plants 



T3SS: type III secretion system; cfu: colony-forming units. 



Functional Genomics Database [28] under accession 
no. D007. 

Flagellin-derived peptides are the major MAMPs causing 
immunity-associated transcriptional changes in tomato 

A first analysis of the genes whose transcript abundance 
was affected by the treatments (Figure 1) revealed two 
main groups based on the overall change compared to 
the mock treatment (>2 fold and P < 0.05). One group 
had a small number of genes (<250; A. tumefaciens, 
DCSOOOAhrcQ-UAfliC) while the other had a much lar- 
ger number (>900; all the remaining treatments). These 
two extreme responses were associated with the presence 
or absence of a perceivable' flagellin (FliC protein) in the 
treatment. FliC in A. tumefaciens is divergent from other 



flagellin proteins perceived by plants, with a non-active 
flg22 region and probably a non-perceivable flgII-28 re- 
gion due to a lack of conserved amino acids in this region 
(Additional file 2: Figure S2). The DC3000AhrcQ-UAfliC 
mutant lacks flagellin but, as with A. tumefaciens, would 
still have other MAMPs. However, our results indicate 
that in tomato other MAMPs play a minor role in gene 
expression changes observed at 6 h. A common feature in 
the second group of treatments is the presence of flagel- 
lin or a portion thereof (flgII-28). In general, for these 
treatments more genes are induced than suppressed, 
and this trend was also observed when protein kinase and 
transcription factor (TF) genes were analyzed separately 
(Figure 1). Even though the titers used for P, fluorescens 
and P, putida were the same, the number of genes affected 



Agrobacterium 
tumefaciens 
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flgII-28 
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2600 



DC3000 
AlircQ-U 
MliC 



DC3000 
AavrPto 
AavrPtoB 



2268 Pseudomonas 
fluorescens 



2484 
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■ I Kinases 

□ D Kinases 

□ I TFs 

□ D TFs 



2531 



Figure 1 Flagellin contains the major MAMPs impacting transcriptional reprogramming following Pseudomonas syringae infiltration 
into tomato leaves. Shown is the total number of genes, protein l<inases (l<inases) and transcription factors (TFs) with increased (I) or decreased 
(D) transcript abundance compared to the mocl<, 6 h after treatment. A >2-fold difference and P< 0.05 were used as cut-offs. The number of 
genes in each category is also shown. See Table 1 for details of the treatments. 
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by P. putida was higher, suggesting a more robust PTI 
response similar to that observed previously in Nicotiana 
benthamiana [29], There are 1,148 protein kinase genes 
predicted in tomato (33% of the 34,727 gene models [30]). 
In our treatments, protein kinase genes represented 
between 7% and 11% of the induced genes, indicating 
an enrichment similar to that observed in Arabidopsis 
upon flg22 treatment [19]. 

We found there was high degree of overlap in the genes 
whose expression changed in response to the flgII-28 pep- 
tide and different bacterial strains with perceivable flagel- 
lin: DCSOOOAavrPtoAavrPtoB (90% overlap), P. fluorescens 
(86%) and P, putida (74%) (Figure 2, see Additional file 3: 
Figure S3 for overlap of suppressed genes). This observa- 
tion, and the fact that we found very little transcriptional 
reprogramming due to DCSOOOAhrcQ-UAfliC infiltration 
indicates that flagellin-derived MAMPs (flg22 and flgII-28) 
have the greatest impact on tomato gene expression changes 
at the 6 h time point. 

The activities of AvrPto and AvrPtoB counteract a 
subset of the transcriptional changes that occur upon 
flagellin perception 

We next performed a clustering analysis of both treatments 
and genes. For this purpose we focused on the 4,150 genes 
that had > 3 -fold difference between the maximum and 
the minimum RPKM observed for any two treatments. 
The genes grouped into two large clusters (Figure 3) each 
with a similar number of genes, with cluster A encom- 
passing genes having decreased transcript abundance 
for flagellin-associated treatments and cluster B having 
genes with the opposite response. Cluster A is enriched in 



gene ontology (GO) terms such as photosynthesis, chloro- 
plast and chlorophyll binding (Figure 3, see Additional 
file 4: Table SI for a complete list of enriched terms). 
In contrast, cluster B is enriched in terms associated with 
defense response, plasma membrane and transmembrane 
receptor protein kinase activity. 

The treatments also grouped into two distinct clusters 
(Figure 3). One cluster included treatments associated with 
flagellin perception (flgII-28, DC3000 AavrPtoAavrPtoB, 
P. fluorescens and P, putida) and the other consisted of 
'flagellin-independenf treatments (mock. A, tumefaciens 
and DCSOOOAhrcQ-UAfliC), Interestingly, DC3000, which 
expresses a perceivable flagellin, clustered with this second 
group. This observation is consistent with the known role 
of AvrPto and AvrPtoB in interfering with FLS2 and 
BAKl, components of flagellin perception [14], and raised 
the possibility that the suppressed genes play particularly 
important roles in the plant defense response. 

The subset of genes induced by flgll-28 but suppressed 
by AvrPto/ AvrPtoB is enriched in gene ontology terms 
associated with the cell periphery, biotic interactions and 
signal transduction 

We focused on the set of genes whose expression is in- 
creased by flgII-28 and is higher in DC3000AavrPto- 
AavrPtoB compared to DC3000 treatment. These are 
genes for which AvrPto and AvrPtoB activity appears to 
counteract the transcriptional increase due to flagellin 
perception. We refer to these as FIRE genes (for flagellin- 
induced, repressed by effectors). From 2,268 genes that are 
increased by flgII-28 treatment (Figure 1), 622 are FIRE 
genes (>2 fold and P < 0.05; see Additional file 5: Table S2). 



flgll-28 



DC3000 MwPt o Mv rPtoB 

' ^ " ^MircQ-U /^fliC 



flgll-28 
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P. putida 
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Total 
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% in 
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with flgII-28 
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DC3000 AhrcQ-U 
AfliC 


223 


67 


30.0 


DC3000 


668 


442 


66.2 


DC3000 AavrPto 
AavrPtoB 


1523 


1364 


89.6 


P. fluorescens 


1706 


1471 


86.2 


P. putida 


2484 


1837 


74.0 



Figure 2 Overlap between genes induced by bacterial strains and flgll-28 treatment. (A) Overlap between Pseudomonos syringoe pv. 
tomato strains and flgll-28; (B) Overlap between soil Pseudomonas and flgll-28. (C) Degree of overlap shown as percentage of genes. 
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1 Category 
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photosynthesis 


4.26E-82 


P 


defense response, incompatible interaction 


2.36E-124 


P 


photosynthesis, light reaction 


1.13E-74 


P 


response to biotic stimulus 


2.69E-124 


C 


plastid part 


1.21E-110 


C 


plasma membrane 


3.64E-61 


C 


chloroplast part 


4.21E-106 


C 


cell periphery 


3.72E-61 


F 


chlorophyll binding 


1.41E-14 


F 


transmembrane receptor protein kinase activity 


6.26E-50 


F 


glucosyltransferase activity 


2.14E-13 


F 


protein serine/threonine kinase activity 


9.49E-49 



I 



■5.00 
■3.33 
-1.67 
0.00 
1.67 
3.33 
5.00 



Figure 3 Clustering analysis reveals that AvrPto and AvrPtoB repress part of the transcriptional reprogramming that is a consequence 
of flagellin perception. Hierarchical clustering was based on genes and treatments. From the 17,597 genes expressed in leaves, those with a 
ratio max/min >3 over the various treatments were included in the analysis. Clusters A and B are enriched in genes decreased and increased by 
flagellin perception, respectively. The color scale shows the deviation from the median for each gene. The top two most enriched GO terms in 
the categories process (P), component (C) and function (F) are listed below each gene cluster. A, Agrobocterium; GO, gene ontology; 
P., Pseudomonas. 



A GO term analysis of these 622 FIRE genes using the 
tomato genome as a reference would not be informative, 
since they derive from a list of genes induced by flgII-28 
and are therefore already enriched in immunity-related 
terms as described for cluster B (Figure 3). For this reason 
we used the 2,268 genes induced by flgII-28 as the refer- 
ence. With this approach we found that the FIRE genes are 
enriched, even within the flgII-28-induced gene subset, for 
immune-related terms and particularly for protein kinases 
and transmembrane receptors (Additional file 6: Table S3). 
A total of 92 protein kinase-encoding genes are among the 
FIRE genes including many that have a known role in the 
immune response: FLS12, Bti9, RIN4J, RIN4.2, Pti4 and 
Lectin RLK VL2 (Additional file 5: Table S2). Furthermore, 
if we had used a less stringent cut-off for defining the FIRE 
gene list (>1.5 fold change and P< 0.05) we would have 
included eight additional immunity-related genes such as 
SOBIRl, SAGlOl and TFT7 (See Additional file 5: Table S2 
for a complete list of genes). 

Given the known importance of protein kinases in the 
immune response, we investigated the families that were 
affected by our treatments (Figure 4). From the 67 pro- 
tein kinase families present in tomato [30], 43 were af- 
fected by flgII-28 treatment. Some families contained 
only genes induced by flgII-28 (for example, ankyrin 
repeat kinases, CTRl/EDRl, MAP2K and MAP3K) and 
others contained only genes with decreased expression 
in response to flgII-28 (for example, LRRK IV and VII 



and RLCK VI). Over 40% of the kinase genes increased 
by flgII-28 (92/214) were found to be FIRE genes. 

Silencing of a FIRE gene encoding a cell wall-associated 
kinase (SIWAK1) compromises pattern-triggered immunity 

To investigate the possible contribution to immunity of 
FIRE genes encoding protein kinases, we selected 33 such 
genes, including at least one member of each kinase 
category, for analysis using virus-induced gene silencing 
(VIGS) in N. henthamiana (Additional file 5: Table S2). 
Silenced plants were subjected to a cell death suppression 
assay (CDSA) [29] to test the involvement of the candidate 
gene in PTI. Briefly this assay consists of an initial infiltra- 
tion of P. fluorescens as an inducer of PTI, then 8 h later 
the infiltration of DC3000 (partially overlapping with the 
first infiltration) as a challenger of this induction. If silen- 
cing of a candidate gene results in a faster breakdown of 
P, fluorescens protection, this is an indication that the gene 
contributes to PTI. VIGS constructs containing a gene 
fragment from FLS2 or an Escherichia co//-derived DNA 
fragment {ECl, see Materials and methods) served as 
positive and negative controls, respectively. From among 
the 33 candidates, we found a cell wall-associated kinase 
(WAK) gene {SlWAKl, Solyc09g014720), whose silen- 
cing compromised PTI-associated cell death suppression 
(Figure 5 A). 5'/ -silenced plants were morpho- 
logically similar to FLS2 and ECl silencing controls 
(Additional file 7: Figure S4). 
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Number of kinases 



10 



10 



15 



20 



25 



30 



35 



Ankyrin Repeat Domain Kinase (7) 
Casein Kinase I Family (16) 
CDC2 Lil<e Kinase Family (20) 
CDPK (38) 
Crinkly 4 Like Kinase (9) 
CRPKl Like Kinase (Types 1 and 2) (11) 
CTRl/EDRl Kinase (13) 
DUF26 Kinase (58) 
GmPK6/AtMRKl Family (25) 
GSK3/Shaggy Like Protein Kinase Family (13) 
IRE/NPH/PI dependent/Se Kinase (30) 
Legume Lectin Domain Kinase (12) 
LRKIO Like Kinase (Type 1) (13) 
LRRK I & Unknown Receptor Kinase I (3) 
LRRKII&X(15) 
LRRK III (11) 
LRRK IV (4) 
LRRK IX (9) 
LRRK VII (6) 
LRRK X (9) 
LRRK XI & XII (62) 
MAP2K (4) 
MAP3K (12) 
MAPK Family (18) 
Other Kinase (7) 
Other Protein Kinase (11) 
Plant External Response Like Kinase (10) 
Putative LRR receptor-like kinase (3) 
Putative protein kinase/Putative RLK (7) 
Receptor-like protein kinase (17) 
RKF3 Like Kinase (2) 
RLCK IV (2) 
RLCK IX (6) 
RLCK VI (9) 
RLCK VII (40) 
RLCK VIII (5) 
S Domain Kinase (Type 1) (17) 
S Domain Kinase (Type 2) (35) 
SNFl Related Protein Kinase (SnRK) (35) 
Unknown Function Kinase (32) 
Wall Associated Kinase (3) 
Wall Associated Kinase-like Kinase (9) 
WNK like kinase - with no lysine kinase (10) 



zz 




□ flgll-28 decreased 0flgll-28 decreased, DC3000 > DC3000 ^avrPto MvrPtoB 

■ flgll-28 increased s flgll-28 increased, DC3000 < DC3000 MvrPto MvrPtoB 

Figure 4 Flgll-28 perception and the activity of AvrPto and AvrPtoB affects gene expression of certain protein l<inase families. Protein 
kinase families are listed on the left with the number of expressed members between parentheses. A >2-fold difference and P<0.05 were used 
as cut-offs. 



To examine the effect of SlWAKl on immunity further, 
we initially performed disease assays using either vacuum 
infiltration or dip inoculation of the N. henthamiana 
pathogen DC3000A/zoj?Q7-7 at 5 x 10^ cfu/mL but we did 
not see any bacterial growth or disease symptom differ- 
ences among £Ci, FLS2- or 5'/ WA/O -silenced plants. As 
an alternative, we first infiltrated plants with P. fluorescens 
at 2 X 10^ cfu/mL to induce PTI and 6 h later infiltrated 
the plants with DC3000A/zoj?Q7-7 at 10^ cfu/mL. Using 
these conditions, we detected reproducible and statisti- 
cally significant {P < 0.01) higher bacterial growth in plants 
silenced for SlWAKl compared with the ECl controls 
(Figure 5B). In ^152- silenced plants bacterial growth was 
intermediate between that of ECl and SlWAKl plants, 
particularly at 3 days after inoculation (dai), but not statis- 
tically different from ECl control plants. Disease symp- 
toms were initially observed at 2 dai, and became more 
apparent at 7 dai (Figure 5C). 5'/ W/l/Ci -silenced plants had 
leaf necrosis similar to FL5'2-silenced plants and both had 
more disease symptoms compared to ECl control plants. 



There are 7 SIWAK and 16 SIWAK-Like (SIWAKL) 
predicted genes in tomato [30]. SlWAKl is a part of a 
clustered four-gene family; two of these are not expressed in 
our conditions (Solyc09g014710 and Solyc09g014730) and 
expression of the other was constitutive (Solyc09g014740). 
A phylogenetic analysis of the WAK and WAKL proteins 
from ArabidopsiSy N, henthamiana and tomato (Figure 6) 
showed that, as expected, those from the solanaceous species 
cluster together. Although some orthologous relationships 
may be evident, the solanaceous WAKs and WAKLs are 
generally distinct from those in Arabidopsis. Interestingly, 
Solyc09g015230 and Solyc09g015240, the tomato WAKs 
that are most closely related to the well-studied Arabidopsis 
WAK genes [31], are not expressed in tomato leaves. 
SlWAKl clusters with several N. benthamiana and tomato 
proteins distinct from the Arabidopsis WAK proteins, 
suggesting the WAKs in these solanaceous species do not 
have a clear orthology to any of the Arabidopsis WAKs. 

We took advantage of the draft genome sequence of 
N. benthamiana [32] and a phylogenetic analysis based 
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Figure 5 VIGS silencing of SIWAKl in Nicotiono benthomiana results in compromised PTI and increased susceptibility to Pst (A) Cell 
death suppression assay. Silenced plants were induced by syringe infiltration of Pseudomonos fluorescens and challenged 8 h later with DC3000 in 
an overlapping area. Cell death symptoms were scored in the overlapping area, on four independently silenced plants with four infiltration sites 
on two different leaves. Results shown are from 7 days after infiltration (dai). Asterisks indicate significant differences based on Dunnett's method, 
using ECl as a control {P < 0.01). (B) Pst growth in leaves. Silenced plants were vacuum-infiltrated with P. fluorescens, inoculated 6 h later with 
DC3000AhopQ/-/ and sampled to measure bacterial populations. Results shown are the average of six plants per construct with the standard 
error. Asterisks indicate significant differences based on Dunnett's method, using ECl as a control (P< 0.01). (C) Disease symptoms of plants 
treated as described in (B) 7 dai. (D) Silencing efficiencies in ECl and SIWAKl silenced plants assessed by qRT-PCR. Silencing efficiency is shown 
as relative expression compared to the ECl control. Data were generated using EFla as reference gene and similar results were obtained using 
PP2A. Data from putative targets: (a) NbSOOOl 1055g0005.1, (b) NbSOOOl 1055g0014.1 and (c) NbSOOOl 1055g0002.1, and a predicted 
non-target (d) NbSOOOl 6938g001 1 .1 are shown. PP2A (e) was used as a cDNA integrity control. Asterisks indicate significant differences using a 
Student's r-test (P< 0.01). dai, days after infiltration; qRT-PCR, quantitative reverse transcription-PCR; Pst, Pseudomonas syringae pv. tomato; PTI, 
pattern-triggered immunity; VIGS, virus-induced gene silencing. 



on nucleotide sequences to identify putative orthologs 
and possible off-targets of the SIWAKl VIGS fragment 
using stringent conditions (that is, genes with >17 nu- 
cleotides of 100% identity). Additionally, we generated 
an RNA-seq dataset from N. henthamiana leaves infiltrated 



with P. fluorescens 55 at 10 cfu/mL, which we used to 
remove genes with low or no detectable transcripts from 
this analysis (< 3 RPKM; Additional file 8: Table S4). This 
approach showed that our VIGS construct was predicted to 
silence three closely related WAK genes in N. henthamiana 
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Figure 6 Phylogenetic analysis of WAK and WAKL amino acid sequences from tomato, Arabidopsis and Nicotiana benthamiana. Green 
font and lines indicate Arabidopsis proteins. Blacl< and red lines indicate N. benthomiono and tomato, respectively. The light gray accession 
numbers for tomato and N. benthomiono are proteins for whose genes few or no RPKM were detected. For all comparisons a >2-fold difference 
and P < 0.05 were used as cut-offs. RPKM, reads per kilobase of exon model per million mapped reads. 



(Additional file 9: Table S5). Quantification of the tran- 
script abundance of these genes by quantitative reverse 
transcription-PCR (qRT-PCR) supported this predic- 
tion as NbS00011055g0005, NbS00011055g0014 and 
NbS00011055g0002 each had decreased transcript abun- 
dance, ranging between 45% and 65%, in plants silenced 
using SlWAKl (Figure 5D). We chose one of the closest 
expressed predicted non-target genes (NbS00016938g0011, 
Additional file 10: Figure S5) and saw no difference in tran- 
script abundance in 5'/ W/l/O -silenced plants (Figure 5D). 
Thus the increased susceptibility to Pst that we observe in 
5/ WA/<7 -silenced plants may be due to either one or a 
combination of three SlWAKl orthologs in N, benthamiana 
but is unlikely to be due to off-target silencing. 

Discussion 

We observed that about 13% of the genes expressed in 
leaves under our conditions were induced 6 h after expos- 
ure to a single flagellin-derived MAMP. We hypothesized 



that Pseudomonas syringae pv. tomato {Pst) has evolved 
to interfere with many of the most important PTI-induced 
genes and so further focused on those genes whose 
expression is suppressed by two pathogen effector proteins, 
AvrPto and AvrPtoB. These effectors target PRR complexes 
as part of their virulence activity and are expected to com- 
promise a large portion of the PTI response emanating from 
subsequent signaling events [14]. Of the 2,268 flgII-28- 
induced genes, 27% were present in this jlagellin-induced, 
repressed by effector {FIRE) class. At least 25 FIRE genes have 
been previously implicated in the plant immune response 
including Bti9 (Solyc02g079600), FLS2 (Solyc02g070890 
and Solyc02g070910), RIN4 (Solyc09g059430 and 
Solyc06g083390), HvCRKl (Solyc03glll540) and LCERK- 
VI2 (Solyc09g005000) (Additional file 5: Table S2). From 
an initial screen of 33 protein kinase-encoding FIRE genes, 
we identified SlWAKl, a gene that, when silenced, com- 
promised the host immune response to Pst infection to a 
similar degree as that observed in FL5'2-silenced plants. 
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Together, the collection of 622 FIRE genes provides a new 
resource for functional genomics screens in solanaceous 
plants and for comparative transcriptomics of PTI and 
ETI responses induced by specific MAMPs and DAMPs 
and diverse pathogens. 

Plants recognize several bacterial MAMPs including 
flg22, EF-Tu, csp22, peptidoglycan and lipopolysaccharides, 
although the importance of each of these in any particu- 
lar plant species varies [4]. Recently, a new flagellin- 
derived MAMP, flgII-28, has been discovered that is 
recognized by tomato and other solanaceous plants but 
not hy Arabidopsis [12]. The possibility that FLS2 might 
be involved in flgII-28 recognition has been excluded 
and its receptor remains unknown [13]. Here we used 
flgII-28 to elicit PTI in tomato leaves and found that it in- 
duces large transcriptional changes similar to host re- 
sponses to other MAMPs. These changes largely overlap 
with those observed when the FliC protein is used to elicit 
PTI as shown by comparisons of wild-type and flagellin- 
deficient bacteria. Together, our RNA-seq experiments in- 
dicate that flagellin-derived MAMPs (flg22 and/or flgll- 
28) are by the far the major bacterial MAMP recognized 
by tomato 6 h post inoculation. We cannot exclude the 
possibility that other MAMPs play an important role at 
other time points, in specific plant tissues or at certain 
stages of pathogen attack. For example, we have reported 
that the Bti9 kinase (the AtCERKl ortholog in tomato) 
likely detects a flagellin-independent MAMP and makes an 
important contribution to Pst resistance, especially in the 
lower leaves of tomato [16]. It also has been discovered that 
in Arabidopsis seedlings the transcriptional response to oli- 
gogalacturonides (OGs) and flg22 greatly overlaps 1 h after 
exposure, but at 3 h the flg22 transcriptional response is 
still high while the response to OGs is diminished, sug- 
gesting that transcriptional reprogramming due to percep- 
tion of these two elicitors changes over time [23] . Further 
experiments with additional MAMPs and at multiple time 
points will be needed to understand better possible tem- 
poral effects on the host response to bacterial elicitors. 

The central importance of flagellin detection in the 
activation of the PTI response in tomato makes it a vul- 
nerable target for bacterial pathogens. Currently, five Pst 
effectors are known to undermine the FLS2-BAK1 path- 
way (AvrPto, AvrPtoB, AvrPphB, Hop All and HopF2) [6]. 
Examination of a 'minimal repertoire' of type III effectors 
needed to confer full virulence to Pst strain DC3000 re- 
vealed that AvrPto and AvrPtoB are early acting' effectors 
whose activity allows for the manipulation by later acting 
effectors of other host defenses such as vesicle trafficking 
[5]. This early activity is due to the fact that AvrPto and 
AvrPtoB each have a domain that binds and inhibits the 
kinase domains of proteins in the FLS2/BAK1 complex 
[14]. Our finding that about 30% of the genes induced 
by flgII-28 are suppressed by AvrPto and/or AvrPtoB 



suggests they may also disrupt the flgII-28 receptor com- 
plex. Each of these effectors also has other domains with 
virulence-promoting activities [15,16,33-35]. Whether these 
domains affect host gene expression is unknown, but 
our present RNA-seq data suggest that transcriptomics 
combined with point mutations in each effector, which 
affect the virulence determinants independently, could 
help elucidate the contributions of each domain. 

Our initial VIGS screen of 33 kinase-encoding FIRE genes 
identified a plant cell wall-associated kinase {SlWAKl) 
that plays an important role in resistance to infection 
by Pst, WAK genes were first identified in Arabidopsis 
where they occur as a clustered five-member gene family 
along with 21 other 'WAK-like kinases' {WAKL) elsewhere 
in the genome [31,36,37]. WAK and most WAKL proteins 
contain a cytoplasmic kinase domain, a transmembrane 
domain and extracellular epidermal growth factor mo- 
tifs [37]. These proteins have been implicated previously 
in host responses to pathogen attack and particularly in 
the perception of pectin-derived fragments [31]. WAKLs 
were found to be highly represented among the kinases 
induced by MAMPs in Arabidopsis although no specific 
gene was shown to contribute to the immune response 
[22]. In fact, a clear demonstration of the role of WAKs in 
host defense by loss-of-fimction analyses has been hin- 
dered by lethality or developmental defects of some WAK 
mutants and possibly redundancy, which is not easily ad- 
dressed by a mutational approach, due to the tight linkage 
of the five WAK genes in Arabidopsis [38] . As a result, the 
best evidence to date that WAKs play a role in immunity 
comes from indirect approaches. For example, overexpres- 
sion of chimeric receptors using domains of WAKl and 
EFR revealed that WAKl is a receptor for OGs and can 
contribute to resistance to Agrobacterium tumefaciens and 
Botrytis cinerea [39]. In another study, a dominant allele 
of WAK2, which was created by fusing it to an artificial 
epitope tag, induced various defense-associated responses 
that were suppressed by a mpk6 mutation [40]. In addition, 
certain WAKLs have been implicated in resistance to fungal 
pathogens [41,42] and VQCQnt[y AtWAKL 10 was reported to 
be expressed coordinately with defense-related genes [43]. 

In our case, silencing reduced the transcript abundance 
of three closely-related WAK genes in N, benthamiana 
by 45% to 65% and resulted in severely compromised 
PTI and enhanced susceptibility to Pst infection. There 
are seven predicted WAK genes in tomato and four of 
them occur in a clustered gene family that includes 
SlWAKl, Unlike genes belonging to the WA/C family in 
Arabidopsis [23], SlWAKl gene expression is induced 
during the response to MAMPs. Extensive sequence 
divergence prevented the assignment of orthologous 
relations between the tomato WAK genes and those in 
Arabidopsis, Consequently, we designated our gene as 
SlWAKl, since it appears to be the first WAK gene to 
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be described in tomato. The previously reported LeWAK 
(Solyc02g090110) [44], is actually a WAKL and is most 
closely related to Arabidopsis WAKL14 (AT2G23450) 
and WAKL21 (AT5G66790) (Figure 6). LeWAK is also 
a FIRE gene and we silenced it along with two other 
FIRE WAK/L genes but observed no effect on PTL 

Based on previous work with Arabidopsis, it is possible 
that SlWAKl acts as a DAMP receptor by playing a role 
in recognition of OGs and subsequent signal transduc- 
tion. OGs could be generated from enzymatic activity 
encoded by either the pathogen or the host [45] . For ex- 
ample, DC3000 has both a pectin lyase (PL) and a poly- 
galacturonase (PG), which could generate OGs during 
the infection process [46]. In addition, tomato has a PG 
gene that was shown to be induced in leaves by damage 
(Solyc08g060970) and, interestingly, this is a FIRE gene. 
Another FIRE gene, Solyc05g005570, is closely related 
to a gene encoding a PG beta subunit that modulates 
PG activity (Solyc05g005560). Inhibitors of PG activity 
have been implicated in contributing to the formation of 
pectin fragments that function as DAMPs [3,45]. A gene 
encoding one of these inhibitors is induced by flagellin, but 
is not a FIRE gene (Solyc07g065090). Collectively, our 
observations support a model in which tomato recognition 
of MAMPs leads to increased expression of SlWAKl and 
an increase in SlWAKl receptors at the plasma membrane. 
Subsequently, OGs are generated by pathogen- and/or 
host-encoded enzymes and are recognized as DAMPs 
by SlWAKl, triggering a sustained immune response. 
Further experiments are needed to test this model, such 
as assaying whether SlWAKl participates in OG percep- 
tion, developing and testing DC3000 mutants that lack the 
PL and/or PG genes, and examining the effect of silencing 
the host PG and PG beta subunit genes. 

Forward and reverse genetic screens for genes that act 
downstream of MAMP recognition have identified rela- 
tively few genes [2,29] . Instead, downstream components 
have mostly been discovered by biochemical (BAKl) or 
transcriptomic (BIKl) approaches [47,48]. This apparent 
shortcoming of genetic screens may be due to functional 
redundancy of specific genes or of PTI-induced pathways. 
It is also possible that, as with BAKl, many PTI-related 
genes play an important role in development or other 
fundamental processes and these genes might be missed 
by forward genetic screens. Here, we excluded from our 
VIGS screen those genes that have already been shown 
to be involved in immunity (25 out of 622). The fact that 
just one gene in the first 33 tested genes had a compro- 
mised immune response might indicate that relatively 
few PTI-induced genes play a critical role in immunity. 
Alternatively, this low frequency could be due to the use 
of heterologous fragments for VIGS silencing, limitations 
of the one PTI assay we used or redundancy in gene 
function. Consistent with the last possibility, multiple 



calcium-dependent protein kinases had to be knocked 
out to observe their involvement in immunity [49]. An 
initial challenge for such combinatorial genetic approaches 
to functional testing of redundant factors is the identifica- 
tion of a set of candidate genes that is small enough to be 
manageable and sufficiently documented to justify the 
effort. The FIRE genes we have identified here meet these 
criteria. Furthermore, our transcriptomic and phenotypic 
analyses of FIRE genes suggest that tomato defense against 
Pst is dominated by perception of a single MAMP and 
one or more DAMPs. Importantly, the FIRE gene set pro- 
vides a new tool for unraveling the redundancies in down- 
stream kinase signaling that are predicted to confer 
system robustness in the face of pathogen effector attack. 

Conclusions 

We performed an extensive RNA-seq analysis and 
VIGS screen to identify novel genes contributing to the 
plant immune response. The host response to flgII-28 
and various bacterial strains indicated that flagellin is by far 
the major MAMP recognized by tomato. A Pseudomonas 
syringae mutant lacking two early acting effectors (AvrPto 
and AvrPtoB) allowed the identification of a defined set 
of FIRE genes. Among these FIRE genes, 25 have been 
previously implicated in plant defense highlighting the 
relevance of this pathogen-defined' set of genes. A wall- 
associated kinase was identified, which, when sflenced, 
compromised the plant immune response to bacterial 
infection. The FIRE genes provide a unique resource for 
dissection of the plant immune system. 

Materials and methods 

Bacterial strains and growth conditions 

Bacterial strains used in this work are listed in Additional 
file 11: Table S6. Pseudomonas strains were grown on 
Kings B medium at 30°C. Agrobacterium tumefaciens and 
E, coli were grown on Luria-Bertani (LB) medium at 30°C 
and 37°C, respectively. Antibiotics used were: ampicillin 
(100 (ig/mL), kanamycin (50 (ig/mL), rifampicin (10 (ig/mL) 
and spectinomycin (50 (ig/mL). 

Plant material and treatments 

Plants were grown at 75% humidity with 16 h light 
(day 24°C, night 22°C). Four leaflets of the third true 
leaf of 4-week-old Rio Grande prfS tomato plants 
were syringe-infiltrated and kept at 24°C and 75% hu- 
midity during the experiment. Infiltrations consisted of 
1 (iM flgII-28 (ESTNILQRMRELAVQSRNDSNSSTDRDA, 
EZBiolab; 90% purity) or bacterial strains (Table 1). 
Nicotiana benthamiana leaves were vacuum-infiltrated 
with Pseudomonas fluorescens 55 {P, fluorescens) at 
10^ cfu/mL and sampled 6 h later. Mock-treated and un- 
treated tissues were also collected. Three biological replicates 
(successive weekly treatments) were performed for each 
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treatment. Tissue was collected at 6 h after infiltration, im- 
mediately frozen in liquid N2 and stored at -80°C until used. 

RNA sequencing analysis 

Total RNA was isolated using TRIzol reagent (Life 
Technologies, Grand Island, NY, USA) and libraries for 
sequencing were constructed as described [50]. Barcoded 
libraries were multiplexed by eight in each lane and 
sequenced on an Illumina Hiseq 2000 system using the 
single-end mode. The length of the reads was 45 or 50 bp. 
Sequence reads have been deposited in the National Center 
for Biotechnology Information sequence read archive under 
accession number SRA096750. RNA-seq reads were first 
aligned to a ribosomal RNA sequence database [51] using 
Bowtie [52] and the aligned read sequences were removed. 
The remaining reads were aligned to the tomato genome 
sequence (version 2.40) using TopHat [53]. Detailed infor- 
mation on the quality of reads in each replicate is provided 
in Additional file 12: Table S7. FoUowdng alignment, for each 
gene model, the count of mapped reads from each sample 
was derived and normalized to RPKM (reads per kilobase 
of exon model per million mapped reads). Differentially 
expressed genes were identified using the DESeq 1.8.3 pack- 
age [54] with the raw count data. Raw P values were cor- 
rected for multiple testing using the false discovery rate [55]. 

For 20% of the genes, we detected 0 RPKM indicating 
these are probably not expressed (or expressed at very 
low levels) in overall leaf tissue in our experimental 
conditions. In a eukaryotic cell with 200,000 mRNA 
molecules, 5 RPKM represents one detected mRNA 
per cell. Leaf tissue is a mixture of different cell popula- 
tions, and the RPKM of a highly expressed but cell-specific 
gene could be masked by the scarcity of the cell type. 
Hence, there is no consensus on how to select and remove 
from the analysis poorly expressed genes, although cut-offs 
ranging from 0.3 to 10 RPKM are found in the literature 
[56,57]. In our case, to set a meaningful RPKM level, we 
examined the RPKM of a subset of genes known to play 
a role in plant immunity, including FLS2 (AT5G046330; 
Solyc02g070890 and Solyc02g070910), Rm4 (AT3G25070, 
Solyc06g083390 and Solyc09g059430), SISERK3A 
(AT4G33430, Solycl0g047140) and Bti9 (Solyc02g079600). 
Of these genes, FLS12 (Solyc02g070910) had the lowest 
values ranging from 0.12 to 6.45 RPKM with 3 RPKM in 
the P. fluorescens treatment (Additional file 1: Figure SI). 
Based on this information we included in our analysis genes 
with >3 RPKM in at least one treatment (approximately 39 
reads mapped per kilobase of transcript for an average of 
13 million mapped reads; Additional file 12: Table S7). 

Virus-induced gene silencing 

For silencing, inserts of 300 to 400 bp were chosen based 
on a BLAST analysis of the N, benthamiana genome 
sequence [58]. Primers were designed using Primer3 [59]. 



PGR amplification was performed using cDNA from Rio 
Grande prj3 tomato leaves treated with flgII-28 (1 (iM for 
6 h; EZ Biolabs, Carmel, IN, USA) or N, benthamiana 
leaves treated with flg22 (1 \iM for 6 h, GenScript, 
Piscataway, NJ, USA). PGR products were cloned into 
pCR8/GW/TOPO (Life Technologies) and recombined 
into a Gateway-compatible TRV2 vector [60]. After 
sequence confirmation, constructs were transformed into 
A, tumefaciens GV2260. VIGS in N, benthamiana plants 
was performed as described [61]. An FLS2 fragment and an 
E, coli gene-based fragment (£C7, which contains a 56% GC 
content and not a single >17 bp-long 100% identical stretch 
in N, benthamiana) were used as positive and negative 
controls. The DNA sequence of the fragments used were: 

>EC7 (Escherichia coli fragment) 

CGGCGTGATTGCGCAAAGCTATCATCAGTCTGAG 

AAATCGGCCTCCGAGTTCGATGCCATTGTTGCGC 

AAACGGAGCAGTTCCTTGCCGACAATGGTCGTCG 

CCCGCGCATTCTGATCGCTAAGATGGGCCAGGAT 

GGACACGATCGCGGCGCGAAAGTGATCGCCAGC 

GCCTATTCCGATCTCGGTTTCGACGTAGATTTAA 

GCCCGATGTTCTCTACACCTGAAGAGATCGCCCG 

CCTGGCCGTAGAAAACGACGTTCACGTAGTGGGC 

GCATCCTCACTGGCTGCCGGTCATAAAACGCTGA 

TCCCGGAACTGGTCGAAGCGCTGAAAAAATGGG 

GACGCGAAGATATCTGCGTGGTCGCGGGTGGCG 

TCATTCCGCCGCAGGATTACGCCTTCCTGCAAGA 

GCGCGGCGTGGCGGCGATTTATGGTCCAGGTAC 

ACCTATGCTCGACAGTGTGCGCGACGTACTGAA 

TCTGATAAGCCAGCATCATGATTAATGAAGCCACG 

CTGGCAGAAAGTATTCG 

>Tomato FLS2.1 (Solyc02g070890) 

AAAGTGTACCGCAGCACTGAGCCTCCAGAGATTT 

TATCAAAAGGATTTGGAACATGCTACCAATAATT 

TCCGTCCGGAAAACATTATTGGAGCCAGCAGTTT 

AAGTACTGTGTACAAAGGAACACTGGAAGATGG 

GAAGATTGTAGCAGTTAAGAAGCTGAATCACCAG 

TTCTCAGCAGAATCTGGTAAATGTTTTGATAGGG 

AAGTCAAGACTCTGAGCCAACTCAGACACAGGA 

ACCTAGTTAAGGTGCTAGGTTACGCTTGGGAAAG 

CAAGAAGCTAAGGGCTTTAGTTTTAGAATACATG 

GAGAATGGGAACTTGGACAACATGATTTATGGTC 

AAGTAGAGGATGACTGGACGTTGTCCAACAGGAT 

TGATATTTTAGTTTCAGTTGCAAGTGGACTATCA 

TACCTGCATTCAGGCTATGATTTTCCAATAGTGC 

ACTGTGACATGAAGCCTTCAAACATTCTTCTGGA 

CAAAAATATGGAAGCACATGTGAGTGACTTTGGG 

ACGGCTAGGATGTTGGGTATTCA 

>SIWAK1 (Solyc09g014720) 

AGGCTACAAACAACTATGCCAGTGATAGAATTCT 
TGGTCGTGGTGGAAATGGAATTGTCTACAAAGGC 
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ATTCTATCTGATAATCGCATAGTTGCTATTAAGAA 

ATCTAAGTTTATGGACGAGGAACAGGTTGAACAG 

TTCATTAACGAGGTACTTATTCTTACTCAAGTCA 

ACCATAGAAATGTTGTGAGACTCTTCGGATGTTG 

TTTGGAAGCCGAAGTTCCTTTACTTGTCTATGAA 

TACATTTCTCATGGAACTCTTTACGAGCATATCCA 

CAATCGAAATGGAGCACCTTGGTTATCTTGGGAA 

AATCGGCTAAGAGTTGCTAGTGAGACAGCAAGT 

GCACTTGCTTACCTTCATTCATCCGCGCAAATGC 

CTATAATTCATAGAGATGTCAAGTCTGCCAATTTA 

TTGTTGGA 

Cell death suppression assay (pattern-triggered 
immunity assay) 

VIGS plants were tested for compromised PTI as described 
[29], with minor modifications. Leaves were induced by syr- 
inge infiltration with P. fluorescens 55 at ODeoo = 0-5 and 
8 h later challenged, in a partially overlapping region, with 
DC3000 at ODeoo^O.Ol. The appearance of cell death 
(CD) in the overlapping area was scored at day 7, using the 
categories: PTI breakdown (>50% of overlapping area with 
CD), and no PTI breakdown (<50%). In each experiment, 
four plants per construct were used with eight overlap- 
ping infiltrations each. Each experiment was repeated at 
least three times. Significant differences are based on 
Dunnetts method, using TRV:EC1 plants as controls 
(P<0.01). 

Bacterial growth assays 

Seven-week-old silenced plants were vacuum-infiltrated 
with a suspension of P. fluorescens 55 (2 x 10^ cfu/mL) 
in 10 mM MgCl2 and 0.002% Silwet L-77 (Comptom Co). 
Plants were kept at 24°C and 75% humidity for 6 h and then 
vacuum- infiltrated with DC3000A/zopQ7-7 (10^ cfu/mL) 
suspended in the same solution. To measure bacterial 
populations, three 0.43 cm^ disks were taken from the 
oldest expanding leaves and processed twice in a Tissue 
Lyser (Qiagen) for 30 sec at 25/sec frequency with 
0.25 mL of 10 mM MgCl2. The volume was adjusted to 
1 mL and serial dilutions were plated on solid LB 
medium with antibiotics. In each experiment, six bio- 
logical replicates per construct were used. Significant 
differences are based on Dunnett s method, using TRV: 
ECl plants as controls {P < 0.01). 

Gene and treatment clustering and gene ontology 
term analysis 

Genes with an expression value of >3 RPKM in at least 
one treatment and a maximum/minimum ratio >3 across 
all treatments were used for clustering. Hierarchical 
clustering was performed using Cluster 3.0 [62] with 
the average linkage clustering method and the Pearson 
correlation as the measure of similarity. GO term 



enrichment analysis was based on the GO::TermFin- 
der module [63]. 

Phylogenetic analysis and degree of overlap between 
bacterial strains and flgll-28 treatment 

A substitution model was evaluated using JModelTest 
[64,65] for nucleotide sequences and ProtTest [66] for pro- 
teins. PhyML through SeaView [67] was used to perform 
the analysis with the default parameters for the GTR model 
(nucleotides) and the JTT model (proteins). One hundred 
bootstraps were used for each analysis. A tree figure was 
created using FigTree [68]. Venny [69] was used to generate 
Venn diagrams to study the degree of overlap between genes 
affected by flgII-28 and different bacterial strains treatments. 

Identification of FIRE immunity-associated genes 

A bibliographic search associated with the FIRE genes was 
performed through: (1) a SQL search of the Sol Genomic 
Network database [70] using the Solyc identifiers in the 
Phenome and Chado schema that contains the manually 
curated loci and (2) a Swissprot database search of the 
Arabidopsis homologous genes (identified using the best 
match of a Selfblast search). The Swissprot database was 
mined using a Perl script (available upon request). 

Virus-induced gene silencing efficiency 

Six N. henthamiana silenced plants per construct were 
vacuum-infiltrated with P. fluorescens 55 as described 
above and samples taken at 5, 6, 7 and 8 h later. Tissues 
from the same plant were pooled and the resulting 
RNA used for quantitative reverse transcription-PCR 
(qRT-PCR). Total RNA was isolated using the Plant 
RNA isolation reagent (Life Technologies, Grand Island, 
NY, USA) according to the manufacturer s instructions. 
Then 6 (ig of total RNA was processed with the TURBO 
DNA-free kit (Life Technologies) for 60 min at 37°C. After 
DNase treatment, 2 (ig RNA was used to prepare cDNA 
using the Superscript III First-Strand Synthesis System (Life 
Technologies) with oligo(dT)20. qRT-PCR was performed 
as described previously [71]. The sequence of the primers 
used for each gene analyzed were: NbS00011055g0005.1 F 
5 '-GAAATATCCCACGGTGATCC-3 NbS00011055g 
0005.1R 5'-GAAGAATGACTGCGGTTAGG-3'; NbSOOO 
11055g0014.1 F 5'-CCCGTTACTCAATACGTTCTT-3', 
NbS00011055g0014.1R 5 '-ATTGGGCGGTGGTTAATG-3 
NbS00011055g0002.1 F 5'-AACAATATCCCACGGTGAC- 
3' and NbS00011055g0002.1R 5'-TTAAAGGAAGACGC 
GAAGG-3'. Cycling conditions during qRT-PCR were 
50°C for 2 min, 95°C for 10 min, and 40 cycles of 95°C for 
30 s, 55°C for 30 s and 72°C for 30 s. Data were normal- 
ized to the NbPP2a and NbEFla genes [72]. The signifi- 
cance of the expression data was analyzed using a pairwise 
Students ^-test {P < 0.01). 
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